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Optimized- Variational Quark Mass Expansion 
and Dynamical Chiral Symmetry Breakdown 



Jean-Loi'c Kneur * 
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A recently proposed variational mass expansion approach to dynamical chiral symmetry breakdown is reviewed. 
We briefly explain how a specific integral ansatz over the analytically continued, arbitrary Lagrangian mass 
parameter, resums the usual variational mass expansion ("delta-expansion" .) The construction is then generalized 
to obtain non-perturbative expressions for the order parameters of the SU{nf)i^ X SU{n[)R breakdown (nf = 
2 or 3), in QCD. Emphasis is put on some general aspects as well as possible limitations of this approach. 



1 Introduction 

To calculate the low-energy properties of the QCD 
spectrum from "first principle" is clearly a desir- 
able task, but obviously out of the reach of or- 
dinary perturbation theory, due to the infrared 
growth of the coupling and the non-perturbativ£ 
relevant dynamics. Chiral perturbation theory EJ 
gives a consistent low-energy effective description, 
but in terms of a set of parameters to be fixed from 
the data, whose precise connection with the basic 
QCD coupling and quark mass parameters is far 
from being resolved at present. Of particular rele- 
vance to the low-energy QCD dynamics are the or- 
der parameters of the SU(nf)L x SU{nf)fi {nf = 
2,3) chiral symmetry breakdown (CSB), such as 
the (qq) condensate or the pion decay constant i<V , 
typically. Partially related to the previous aim, it 
has however been known for a long time that, at 
least in simplified ficld-thcorctic models, definite 
information on specific non-perturbative quanti- 
ties may be inferred from particular resummation 
propertiesci and/or appropriately modified pertur- 
bation seriescL 

Recently Qu we have investigated a new ap- 
proach to explore how far the basic QCD La- 
grangian can provide non-zero quark condensates, 
pion decay constant, as well as dynamical quark 
masses, in the limit of vanishing Lagrangian (cur- 
rent) quark masses. The starting point is very 
similar to the idea developed long ago in refs. u, 
where the convergence of ordinary perturbation 
was shown to be improvable by a variational pro- 
cedure in which the separation of the action into 
"free" and "interaction" parts is made to depend 
on some set of auxiliary parameters ^. 

" In one-dimensional field theories this optimized per- 



An essential novelty ialhat our construction can 
reach infinite order uu of such a variational- 
perturbative expansion, therefore presumably op- 
timal, provided it converges. It also gives a consis- 
tent treatment of the renormalization, reconciling 
the variational expansion with the inherent infini- 
ties of most field theory of dimension higher than 
1. 

2 Mass Gap in the Gross-Neveu model 

We shall first illustrate the basic ineijedients of our 
construction with a determination Lru of the mass 
gap in the 0{2N) Gross-Neveu (GN) model. 

From renormalization group (RG) resumma- 
tion properties, one can infer the following form of 
a bare, RG- invariant resummed mass (D = 2 + e): 

mF{mo)^mo{l~4,TTbogQTmpy^ , (1) 

where mo, go are the bare mass and coupling, T = 
r[-e/2]/(47r)i+^/2, and 6o, 7o are the one-loop 
RG-coefficients of the running coupling and mass, 
respectively (we use a normalization such that 

Pig) = -bo9^ - hg^ , imig) = log^ + 715" + 

• • •; where in the 0(2AD model 60 = (A^- l)/(27r), 
70 = (AT _ l/2)/7r, etcH.) As easily shown, eq. (||) 
gives, after renormalization, the exact mass gap 
result in the A ^ 00, to ^ limit: mp — A-jg^, 
where A-p^- = /x e"-'^/''^''"^^^'^)^ is the basic scale 
in the MS scheme. Now, in the more compli- 
cated case of arbitrary A^, according to the above- 
mentioned variational principle, we expect to ob- 
tain a series of approximants to the mass gap, 

turbation theory ("delta-expansion") has been shown E to 
lead to a rigorously convergent series of approximations, 
even in strong coupling cases. 
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by optimizing with respect to the arbitrary La- 
grangian mass, at successive orders of the vari- 
ational a;-series, formally defined by substituting 
everywhere in the bare Lagrangian: 



mo ^ Too (1 - x); 



9l 



X, 



(2) 



and in particular for example in expression (0). 
However, in order to get a finite and non-trivial 
result (i.e. mp ^ Q) one has to resum the resulting 
a;-series by using a specific contour integration^ of 
(|l|) over the (analytically continued) x parameter. 
Moreover, the one- loop RG-resummed expression 
(|l|) , which includes only the leading dependence in 
TO, has to be generalized to include both higher or- 
der RG-dependence upon to, and non-logarithmic 
purely perturbative corrections as well, in order to 
get a more realistic mass gap in the arbitrary N 
case. In this way we end up with the ansatz: 



(to ) 



A' 



MS 



2-^to" 
2m 



dv- 



FA{v)[C + F{v)]B 
Ml M2 



(3) 



F{v) F^{v) 

where the contour is around the negative real axis, 

F{v) = \ii[m"v]-A \nF-{B-C) ln[C+F], (4) 

with A - 7i/(26i), B = 7o/(26o) - 7i/(2&i), 
C = &i/(26§); A-j^jg is the (RG-invariant) scale at 
two-loop order; for convenience we introduced the 
scale-invariant, arbitrary (dimensionless) "mass" 
parameter, redefined in terms of to as 



_ / TO(/i) 



A 



MS 



2^ [2&05 



2i-7o/(26o) 



1 + ^.9^ 



n B 



(5) 

As mentioned, F{1) in the integrand of (|^) by con- 
struction resums the leading and next-to-leadiiig 
logarithmic dependence in to(/2) to all orders 
The two-loop non-logarithmic perturbative coeffi- 
cients M.1, A^2, given explicitly in ref. connect 
the pole mass with the running mass m{p,). 

We have shown □ that expression (||) does re- 
sum the x-series [| generated from the substitution 
(^. Moreover it is possible to choose a renormal- 
ization scheme (RS) such that bi = 'ji — for 
i > 2. In that sense, eq. (|^) resums the full RG 

''v in is related to the original expansion parameter 
xasx = l — v/q, q being the order of the i-expansion. 



dependence in x and to". In contrast, the purely 
perturbative (non-logarithmic) information, con- 
tained in Ail, Ai2, is limited to the two- loop or- 
der. This is where the variational principle and 
optimization play their role, whereby we hope to 
obtain a sensible approximation to the exact mass 
gap. Accordingly, we may now look for extrema 
of expression (|^) with respect to to", using stan- 
dard contour integration techniques. Observe in 
fact that, were we in a simplified theory where 
Ali=A^2 = -- ' = (as incidentally is the case in 
the large- limit), (^ would have a very simple 
behaviour near its optimum (at to" — > 0), giving 
a simple pole with residue M2 = {2C)~^ ^ms- 
Now, in the arbitrary N case where A4i, Ai2,--- 
cannot be neglected, one can construct a set of ap- 
proximants of (^) in the to" limit, with som£j 
variants of Pade approximants (PA). The resultscl 
of two different optimal fifth-order PA are com- 
pared in the table below with exact results for the 
mass gap of the 0{2N) model (obtained alterna- 
tively from the Bethe ansatz in ref. E3). As illus- 
trated, an error of 0(5%) or less, depending on N , 
can be obtained |. 





Exact 


Pade 1 


Pade 2 


N 


tof/A^ 


(error) 


(error) 


2 


1.8604 


1.758 (5.5%) 


1.875 (0.8%) 


3 


1.4819 


1.475 (0.5%) 


1.486 (0.3%) 


5 


1.2367 


1.284 (3.7%) 


1.265 (2.3%) 



3 Dynamical quark masses 

Most of the previous construction (for arbitrary 
N) may be formally extended to the QCD case. 
In particular, expression (^), with the appropriate 
change to QCD values of the hi and 7^ RG,|Coeffi- 
cients, provides a dynamical mass ansatz QU as a 
function of the quark flavours n/ (the latter en- 
tering both A-jgg and the bi, 7^ expressions). It 
is important to note that expression (|3|) for arbi- 
trary N in the GN model uses exactly the same 
amount of perturbative and RG information as is 
at our disposal at present for a QCD quark mass: 
namely, the exact two-loop_B,G-resummed plus 
purely perturbative A^i, A42c^ dependence. Since 



'^The results given in the table were obtained by further 
optimizing with respect to an arbitrary scale parameter, 
introduced from /i — > a /i. Results from PA of order lower 
than the (optimal) order 5 give a larger error. 
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our construction only relics on RG-properties (and 
analytic continuation techniques), passing from 2 
to 4 dimensions is not expected to cause major 
changes, at least naively 0. 

Actually, one complication does occur: as a 
more careful examination of relation (Q) indicates, 
there are extra branch cuts in the v plane, with 
Re[wcnt] > for the relevant case of n/ = 2 or 3 
in QCD. This prevents using the expansion near 
the origin, which would lead to ambiguities of 
C'(exp(+const./m")) for m" 0. The origin of 
those singularities is rather similar to the renor- 
malon onesB, since they appear in a resummed ex- 
pression relating a "reference" scale Mdyn — ^jjs 
to an infrared scale m" ~ 0. However, in the 
present construction it is possible □ to move those 
extra cuts to the safe location Re[Uj,„(] < 0, by 
noting that the actual position of those cuts de- 
pend on the RS, via 71. Performing thus a second- 
order perturbative RS change in m{^), g{fJ.), which 
changes 71 — > Ji, and looking for a flat opti- 
mum (plateau) of an appropriately constructed 
PA, with respect to the remnant RS arbitrariness'!, 
we obtain: 

<pf'K'-0):^2.97A^(2) (6) 
for Hf — 2, and a similar result for Uf — 3. 



4 Order parameters of CSB: Fj, and (qq) 

The above dynamical quark mass, although it 
has some meaning as regards spontaneous CSB in 
QCD, hardly has a direct physical intepretation, 
e.g. as a pole in the S-matrix, due to confinement. 
In other words, it is not a properly defined order 
parameter. It is however possible to generalize 
the mass ansatz (^) to obtain a determination of 
the ratios F.^/Ajjg and (gg)(^)/A|— . The latter 
gauge-invariant quantities are unambiguous order 
parameters, i.e. 7^ or (qq) ^ indicate spon- 



From that point of view, there are no differences be- 
tween QCD and the GN model (for arbitrary N), as com- 
pared e.g. with other two-dimensional (or a fortiori, one- 
dimensional) theories where particularly simple features 
are sometimes due to their super-renormalizable properties. 
Note also that a mass gap in the (pfflAf) GN model does not 
contradict the Coleman theorem til on the non-breaking of 
the continuous chiral symmetry in two dimensions. This in 
a sense illustrates that the applicability of our method is 
more generalithan the physics of CSB. 
''See ref.El for details. 



taneous CSH. The appropriate generalization Q of 
(|) for isB 



i?2 

A^ 

MS 



(2&o) 



-2C 



ll\2 



2i7r 



dv 



5^ 1 



p2A-l'^C + F] 2S 

in terms of the very same F{v) defined in eq. (^) 
(therefore leading to the same extra cut locations 
as in the mass case), and where 5t^, a-^r and 
are fixed by matching the perturbative MS ex- 
pansion, known to 3- loop orderllj'cj. A numerical 
optimization with respect to the RS-dependence 
of an appropriate PA, along the same line as the 
mass case, gives e.g for 71/ = 2: 



Fi'^f (m" 



.,optv- 0) ^ 0.55 A^(2) . (8) 

Concerning {qq), an ansatz similar to (Q) can be 
derived (with coefficients 5, a, (3 specific to {qq) 
and obvious changes in the m", F and v powers), 
but for the RG- invariant combination m{qq) only, 
since our construction only makes sense for RG- 
invariant quantities. As it turns out, the limit 
TO ^ simply gives m{qq) — > 0. To extract 
an estimate of the (scale-dependent) condensate 
{qq){fi) irL-our framework is only possible by in- 
troducing B an explicit symmetry-breaking quark 
mass niexp (independent from to), and expanding 
the m{qq) ansatz to first order in rriexp- This gives 
e.g. for Uf — 2: 

(q-q)(A2 = lGeV)c^0.52A^(2). (9) 

Comparing the results (^), ^ and (^ gives a 
fairly small value of the quark condensate (and a 
fairly high value of the dynamical mass) , as, com- 
pared to other non-perturbative methods 113. Al- 
though small valueSpOf {qq) are not experimentally 
excluded at present it is also clear that our rel- 
atively crude approximation deserves more refine- 
ments for more realistic QCD predictions. 

5 Summary and conclusions 

The variational expansion in arbitrary to, devel- 
oped in the GN model eI, can thus be formally ex- 
tended to the QCD case. This gives non-trivial 

non-trivial point in this generalization is the nec- 
essary additional (additive) renormalization of the axial 
vector-axial vector two-point correlator, whose first-order 
exfransion term with respect to external momentum defines 
FttU. For a fixed RS, this unambiguously leads to (M. 
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relationships between ^jjg and the dynamical 
masses and order parameters, and (qq). 

Let us conclude with some remarks on the 
possible generalizations (and limitations) of this 
approach. In principle, it may be applied to any 
(renormalizable or super-renormalizable) field the- 
ory, where the massless fermion (or scalar) limit 
is of interest, and most probably the basic idea 
could be applied to parameters other than mass. 
Now, what is at any rate limitative, is the rela- 
tively poor knowledge of the purely perturbative 
part of the expansion (only known to two-loop or- 
der in most realistic field theories). Our ansatzes 

(0) would be exact if there were no such cor- 
rections, i.e. if the dependence on m was entirely 
dictated from RG properties. Therefore, our fi- 
nal numerical results crucially depend on the opti- 
mization |[ Apart from a few models where the se- 
ries is knospito large orders (as in the anharmonic 
oscillator or in the GN model for N — > 00), 
we can hardly compare successive orders of the 
variational-optimized expansion to estimate, even 
qualitatively, its convergence properties. Invoking 
the "principle of minimal sensitivity" 0, although 
physically motivated, in a sense artificially forces 
the series to converge, with no guarantees that it 
is toward the right result. 
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